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^ \ Abstract 

A global existence result is established for a free boundary prob- 
lem of planar magnetohydrodynamic fluid flows with radiation and 
^ , large initial data. Particularly, it is novelty to embrace the constant 

■ transport coefficient. As a by-product, the free boundary is shown to 

. expand outward at an algebraic rate from above in time. 

o 

1 Introduction 

In astrophysics, stars may be viewed as compressible fluid flows formulated 
by the Navier-Stokes equations, which is now expressed by the conservation 
^ '. of mass, balance of momentum and energy. However, their dynamics are of- 

^ \ ten influenced by magnetic fields and high temperature radiation. When the 

radiation is taken into account, the new balance law for the intensity of radia- 
tion should be added to complete the hydrodynamic system. More precisely, 
the material flow needs a relativistic treatment since the photons are massless 
particles traveling at the speed of light. In other words, it raises more difficul- 
ties in mathematics and physics. Fortunately, under some plausibly physical 
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hypotheses and asymptotic analysis, especially for the equilibrium diffusion 
model, the radiative transfer equation can be replaced by the usual hydrody- 
namic model equations with the pressure, the internal energy and the heat 
conductivity added by the special radiative components, see for instance [15]. 
Furthermore, when the magnetic fields are considered, the motions of con- 
ducting fluids may induce electric fields. Thus, the complex interaction of 
magnetic and electric fields significantly affects the hydrodynamic motion of 
fluid flows, which is described by Maxwell's equation. Finally, the system of 
magnetohydrodynamics is the following in Eulerian coordinates. 

dtp + div(pu) = 0, (1.1a) 

dtipu) + div(pu (S)u) + Vp = divS + (V x H) x H, (1.1b) 

dt{pS) + div{pS'u + pu) + divQ = div(§u + z/H x (V x H)) 

+ div((u X H) X H), (1.1c) 

Ht - V X (u X H) = -V X (z/V X H), divH = 0, (l.ld) 

where p is the density of fluid flows, p is the pressure, u G the velocity, 
H G M'^ the magnetic field. S stands for the total energy expressed by 

£ = p(e + i|uP)+i|HP. £' = p(e + i|„|=). 

which e is the internal energy of flows. The viscous stress tensor 

§ = A'(divu)I + fi{Vu + Vu^), 

where A' and p are the viscosity coefficients of the flows with X' + 2p > 0. 
I is the 3x3 identity matrix. Vu""" is the transpose of the matrix Vu. 
Particularly, we emphasis that the electric field is induced by the velocity u 
and the magnetic field H, 

E = z/VxH-uxH, 

and the induction equation (ll.ldl) is derived from neglecting the displacement 
current in Maxwell's equation. For more detailed physical explanation and 
mathematical deduction, see the appendix of [Ij for reference. 
The pressure p of the gas obeys the equation of state 

P = PG+PR = Rp0 + l0'- (1.2) 
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R is a constant depending on the material properties of the gas, 6 the tem- 
perature of fluid flows and the last term in the above equality denotes the 
radiative pressure with a > being the stefan-Boltzmann constant. Accord- 
ingly, the internal energy is 

e = eG + eR = CJ + ^9^ (1.3) 

with Cv being the heat capacity of the gas at constant volume, and similarly 
the heat flux Q is 

Q = _kV0 = - (/^i + V^, 

where ki is a positive constant and c is the speed of light. The quantity 
1/kp is the mean free path of a photon inside medium, which is related to 6. 
Without loss of generality, we can assume that 

+ e'^)<K = k{p, 9) < + 6''), g > 0, (1.4) 

where ni, H2 are positive constants. 

In the present paper, we primarily study a free boundary problem for 
planar magnetohydrodynamic fluid flows. More precisely, the motion of flows 
is supposed to be in the x-direction and uniform in the transverse direction 
{y,z), i.e., 

p = p{x,t), e = e{x,t), 

u = (m, w)(x,t), W=(m2,M3), 

H= (6i,b)(x,t), h = {b2,h), 

and the corresponding dynamic equations (II. ip are reduced to the following 
in the Eulerian coordinates. 

Pt + {pu)cc = 0, xent:= (a(t), b{t)), t > 0, (1.5a) 

{pu)t + {pu^ +p+^\h\^)^ = (Am^)^, (1.5b) 

(pw)i + {puw - b)^ = (/iw^.)a;, (l-5c) 
ht + {uh - w)^ = (z/b^)a;, (1.5d) 
{pe)t + (peu)^ - {k9^)^ = Xul + /i|w^p + i'\h^\^ - pu^. (1.5e) 



3 



where A = A' + 2/i and 61 = 1. 

To supplement the system f ll.Sp . we impose the following initial condi- 
tions, 

(p,M,6', w,b)(x,0) = (po, Mo, 6^0, Wo, bo) (x). (1.6) 
and the free boundary conditions 

(AM,-p,w,b,^,)(/(t),t) = 0, f{t)=a{t), b{t). (1.7) 

where f{t) denotes the free boundary defined by f'{t) = u{f{t),t). 

The aim of this paper is to establish the well-posedness of the system 
fll.5p - fll.7p for the general setting of the heat conductivity fll.4p . especially 
including the case of constant transport coefficients. 

Let us first review some related works in this line. For instance, Zhang- Xie 
[28] studied the global existence of the equations f ll.Sp for Dirichlet boundary 
problem when g > |, which is extended to g > (2 + V211)/9 by Qin-Hu 
[T5] . In the forthcoming companion paper pT], we also showed that global 
solutions exist for large initial data even if g > 0. On the other hand, Chen- 
Wang [HEl] established the existence of global solutions of real gas for a 
free boundary or Dirichlet problem, respectively, where they assume that the 
pressure and internal energy satisfy the following condition with exponent 
r e [0, 1] and g > 2(r + 1). 

0<pp<Po{l + e'+'), ee>eo{l + en. 

We remark that the restriction on g in fll.4p is only from a more mathe- 
matical point of view. The main difficulties come from the interaction of 
magnetic fields and fiuid velocity, as well as higher nonlinearity of radiative 
term, which prevent solving the initial boundary problem by means of known 
analytic techniques and tools. As a consequence, some essential new ideas 
are proposed to track this open problem. 

From the previous works, we know that one of the key points of the 
problem is to achieve the upper and lower bounds of density and temper- 
ature, which implies that there is no concentration of mass and heat. By 
delicate analysis, we notice that the radiative term helps us get the up- 
per and lower bounds of density based on entropy-type energy estimates. 
Due to the interaction of the dynamic motion of the fiuids and magnetic 
field, a priori estimates of temperature are much more complex. In our 
context, new a priori estimates of 1 16*^1 1 2,1 [0,1] and | |px| |l2[o,i], which are con- 
trolled by \\yVxx\ \Lmo, i]xio,t]) and | |b^x| U2([o,i]x[o,i]), are proposed, see Lemma 
13.91 and Lemma fS.lOt and then a priori estimates of | |wa,3,| |i2([o,i]x[o,t]) and 
I Ibxil |L2([o,i]x[o,i]) are subsequently obtained. With these bounds, all required 
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a priori estimates for 1 16'^| |ii[o,i], | IPxI |l2[o,i] and | |L4([o,i]x[o,t]) are achieved, 
refer to fl3.24l) - fl3.26l) . In the subsequential process, motivated by p!|[TT|[T3]. 
we succeed in getting the upper bound of temperature and the first derivative 
of velocity. 

To our best knowledge, there is a few results for the case of perfect flows, 
namely neglecting the radiative effect, or equivalently a = in fll.2l) . when the 
transport coefficients are positive constants. Among them, Kawashima and 
Okada [12] proved the existence of global smooth solutions to one-dimensional 
motion with small initial data. In addition, Hoff and Tsyganov [8] consid- 
ered uniqueness and continuous dependence on initial data of weak solution 
of the equations of compressible magnetohydrodynamics. Moreover, Fan- 
Jiang-Nakanura [7] showed the global weak solutions of plane magnetohydro- 
dynamic compressible flows converge to a solution of the original equations 
with zero shear viscosity as the shear viscosity goes to zero when q > 1. 
However, it is an open problem for the global existence of perfect flows with 
large initial date, which may be shed light on by our techniques. From our 
analysis, we know that the focus of the problem is on the lower bound of 
density since the lower bound of density in our context is strongly depended 
on the radiative term. 

However, similar obstacle on q is also in existence even if we neglect 
the magnetic fields. For instance, Ducomet-Zlotnik [3] proved the existence 
and asymptotic behavior for ID radiative and reactive gas when q > 2 and 
Umehara-Tani [221 [23] made further extension in this direction for ID or 
spherically symmetric case when 3 < q < 9, which extended to g > by 
authors in [16] . In addition, Wang and Xie [27J showed global existence of 
strong solutions for the Cauchy problem when q > ^ and the reference [1] 
proved global-in-time bounds of solutions and established it global exponen- 
tial decay in the Lebegue and Sobolev spaces when g > 2. 

Indeed, there are also extensively studies on MHD hydrodynamics, which 
is beyond our ability to address exhaustive references, see for instance [5l[6l 
[9l[l0l[l9H2ll[25l[26l[2H] and references cited therein. 

We formulate our problem and state the main results in section 2, and In 
section 3, we give the essential a priori estimates for the existence of global 
solutions. 

2 Preliminaries and Main Result 

Before stating the main result, we first introduce the Lagrangian coordinates 
to translate the free boundary problem fll.5p - fll.7l) to the fixed one for con- 
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vience, which is equivalent under consideration. Let 

\{U)dt t = t. (2.1) 

a{t) 

Then < y < 1 = Jq p{^, t)d^ which is the total mass of fluid flows without 
loss of generality. The system (11. 5p canonically becomes 

vt = Uy, (2.2a) 

ut=(-p-hh\' + ^) , (2.2b) 



2' ' V 



y 



w.= (b+^i^), (2.2c) 

V V Jy 

{vh),= (w+'^) , (2.2d) 



' y 



A \ /i|wy| z/|b. 



where f = 1/p is the specific volume. The corresponding initial data and 
boundary condition are as follows 

(w,M,^,w,b)(?/,0) = (t;o(y),Mo(l/),^o(2/),wo(?/),bo(?/)). (2.3) 

and 

(A^-p,w,b,^^^)|(ci,t) =0, rf = 0,l. (2.4) 

With these preliminaries, we are now ready to state the main result for 
the system (l^ - fl^ . 

Theorem 2.1. Let g > and a G (0,1). Suppose that A,/i, z/ are positive 
constants. In addition, assume that the initial datavo{y), uo{y) ,wo{y) , ho{y), 
9o{y) satisfy 

Co'<voiy), eoiy)<Co, 
for some positive constant Cq and 

{voiy),uoiy),woiy),hoiy),Ooiy)) e C'+'^in) x C'+-{n)<', 

for a G (0, 1). Then there exists a unique solution {v, u, w, b, 6) of the initial 
boundary value problem (12.21) - (12.41) such that 



{v,Vy,vt)eC'''''/\QTY 
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and 

Moreover, the expand rate of interface is 

< b{t)-a{t) < C{l + t), 
where C is a positive constant independent of time t. 

Remark 2.1. Obviously, the case of constant transport coefficients is in- 
cluded when q = 0. In addition, it is also valid for n = Ki + K2— instead of 

The procedure for the existence of global-in-time solutions is classical 
from the standard method by the global a priori estimates of (p, u, 6, w, b). 
Therefore, the main task for us is to establish the global a priori estimates. 

3 Some a priori estimates 

In order to establish the existence of global solutions, we need to deduce some 
a priori estimates for (p, u, 6, w, b), which is essential to extending the local 
solutions by fixed point theorem. In the sequel, the capital C{C(T)) will 
denote generic positive constant depending on the initial data (time T > 0), 
which may be different from line to line. 

3.1 Upper and lower bounds of density 

As usual, the basic energy estimates is as follows. 
Lemma 3.1. We have 

[e + ]^{u' + |w|2 + v\h\^)^ dy < C. (3.1) 
Proof. From fl2.2ap - fl2.2ep . we get from the boundary conditions (12. 4p 
|^'(e + i(,.' + |wP + „|br))<iy = 0. 

which implies (13. ip . □ 

And then, we immediately arrive at the uniformed upper boundedness of 
density for all time. 
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Lemma 3.2. Suppose the hypotheses of Theorem \2.1\ are valid, the 

p{y,t)<c, 

and 



and 



e^dy < C, 



m'^dy < C. 



for (y,t)G(0,l)x(0,T). 
Proof. The equation fl2.2bl) can be rewritten as 

Ut = ^X{liav)t -p - ^|b|^ 

Integrating it over [0,y] x [0,t], we obtain 

Alnf = Alnfo(?/) + / {p + -\h\'^)ds + I {u — uo)dx 



and notice that 



(u — Uo)dx 



•1 rl 

< C I u^dy + C uldy < C, 
'o Jo 



Inw > -C, 



which imphes 
or equivalently 

p{y.t)<c. 

Following from fl3.ip and f l3.2p . we easily get 



1 i-l rl i-l 

e^dy = / pve'^dy <C vO^dy <C edy < C. 
Jo Jo Jo 



Similarly, one has 



1 /■! p1 

2 J„. _ / „„.lul2j„. / / „.lul2. 



Jo Jo 

This ends the proof. 

In the sequel, the expand rate of free boundaries is obtained. 



8 



Lemma 3.3. It holds that 

< / vdy< C{l + t), 
Jo 

or equivalently, by the transformation (12. ip . one has 

< b{t)-a{t) < C{l+t). 
Proof. Integrating [0,y] over ( 12.2b|) . it yields 

\tdx = ^-{p + l-\h\'), 

Q V Z 

which imphes by multiplying v and integrating with respect to y and t 
X J vdy 

1 rt /•! 

2 



A / Vo{y)dy+ I / v[p+-\h\ ] dyds 



^0 ^0 

t pI 




^0 



V I / Utdx ] dyds. 



For the second term on the right hand side, it is deduced that 




v[p+ ^|b|' ) dyds <C [ [ (e + v\h\'^) dyds < C{1 + t). 
'0 Jo \ ^ J Jo Jo 

On the other hand, we get from f l2.2ap by integration 

v{y,t) = Vo{y) + Uy{y,s)ds. 
Jo 

Thus, it leads the third term on the right hand side to 
V / Utdx dyds 




Jo 



(^vo{y) + j Uy{y,s)ds^ udx^dy 

+ / v^dyds — I vo{y) I / uo{x)dx j dy 
Jo Jo Jo \Jo J 

vo{y) {u — uo{x))dx^ dy + J J v^dyds 



u{y,t) u{y,s)dsj dy, 

< 0(1 + t). 
Thus, we utilize (13. ip to finish the proof. 



9 



Lemma 3.4. One has 



where 

U{t) 

and 



U{t)+ fv{T)dT<C{T), (3.7) 
Jo 

[ {C,ie-l-\oge) + R{v-l-\ogv))dy, 
Jo 



v(t)= + ^ + ^ + 

forO<t<T. 

Proof. According to the expression of p and e, we thereby compute 

eg9t + OpgUy + + -^-^ +{-dy) , 

V V V \V )y 

which imphes that 

(c^Xogd ^RXogv^^avQ^dy 



t;^ t;^ t;^ vQ'^ ^ 
Integrating that over (0, 1) x (0,t), it yields 

Vit) + ^ r(r)rfr <C (^1 + vO'^dy 

On the other hand, one has by (13.61) 



ve^dy < v9^dy^ ' (^j^ vdy 
< C{T), 

which ends the proof in conjunction with (13. ip . □ 

Furthermore, we can deduce some a priori estimates for 6 and b for all 
time. 
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Lemma 3.5. We have 



f max^«+^(y, s)ds < C(T), g > 0, 

Jo [0.1] 



and 



(3i 



(3.9) 



forO<t<T. 

Proof. There exists y{t) by the mean value theorem such that 



o{yit),t)= f'edy, 

Jo 



where y{t) G [0,1] for each t G [0,T]. Furthermore, it yields by Holder 
inequality 



0{y,t) 



q + i 



edy 



q+4 



+ 



q + A [y 



q + i 



y{t) 



e{^,t)—-%{^,t)d^ 



<c 1 



< c 



< c 



K2 



1 yQI+i 



—dy 



K 

1 yQq+i 



dyj V{t)-2 

1 

dyVvit)'^ 







< C 



1 + 



ve^dy ) V{t)^ 



< C (l^Vit)^ 



Then taking square on both sides of the above inequality, we can get (13. 8p 
by integrating it over [0,t] with the help of (13. ip and (13. 7p . 

The proof of (13.91) is directly deduced from (13. ip and (13.80 . Indeed, re- 
calling the boundary conditions (12. 4p . we find that 



b|2 = 2 / h-hi.{i,t)di 



< C 



v9 



dy + C ve\h\'^dy, 
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or 



I 



^\\L°^no i])ds < C / / -^dyds + C / max6' I / t>|b| j (is. 

Jo Jo "^^ io [O'l] \Jo 



< C{T). 

This finishes the proof. □ 

With these preliminary processes, the lower bound of density can be in- 
ferred that 

Lemma 3.6. One has 

p{y,t)>C{T). (3.10) 
Proof. Recalling the equality (13. Sp . we find that 

ft I ry 
Xhiv = XhiVo{y) + / (p + -|bp)(ir+ / {u — Uo)dx 



2 JO 

t rt fl 

2 j„ , n I „,2. 



<C{T) + C mixxO'^ds + C I |b| |ioo([o,ii)C/s + C / u'dy 
Jo [O'l] Jo ' Jo 

< C{T), 

where we have used Lemmas 13.11 and 13.51 □ 

3.2 Upper and lower bounds of temperature 

At first, we are in a position to deduce a priori estimates of the derivatives 
of (m, w, b). 

Lemma 3.7. One has 

[ [ {ul + \wy\^ + \hy\^)dyds<C{T). (3.11) 
Jo Jo 

Proof. The equality fl2.2ep can be rewritten as 

„H-p„„ = M + ^lK^ + !l!M+(!i»,) . (3.12) 

V V V \V Jv 
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Then integrating fICT]) over [0, 1] x [0,t], we get from ([XS]) and ([SS]) 




JO 



V V V 

(e - eo)dy + 



npUydyds 



<C + - I I —^dyds + C I I p^dyds 

<C + 



Jo 'J Jo Jo 



-I I ^dyds + C [ maxe'^f [ O^dy] ds 
2 Jo Jo V Jo m \Jo J 



2 _ 

< C(T) + - / / -^dyds, 

2 Jo Jo V 

which ends the proof according to (13. 8p and (13.101) . □ 
Lemma 3.8. The following inequalities hold for b 

f [ \h-h,/dyds<C{T), (3.13) 
Jo Jo 



and 

ft rl 




\h\^dyds < C{T). (3.14) 



^0 

when t G [0, T]. 

Proof. By (13. 4p . one has 

ft ^1 

\hfdyds 

t / rl 




Jo 



< [ maxlbl^ ( [ IbPrfu^ ds 

- Jo m \Jo J 

<C [ [ \hfdyds + C [ [ \h\^\h ■ hy\dyds 
Jo Jo Jo Jo 



^0 ^0 ^0 



<-/ / m^'dyds + C I I ^ — ^^dyds + C{T), 



which imphes 



t rl ft rl |U|2|U 12 

\hfdyds<C / ' ' ' dyds + C{T). (3.15) 




^0 ^0 JO 



V 
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Multiplying 4|bpb on both sides of fl2.2dp . one has 

(v\h\% = + ■4|b|2b-3^;,|b|^ 

which follows from fl3.9p and (13. lip that 

-1 rt /•! IKI2IK 12 



r bP b P 

/ vlhl^'dy+Uu / ' ' dyds 
Jo Jo Jo V 



1 fi i*t /»i 

v\h\^{y,{i)dy — 12 I / \h\^w -hydyds — ?> I / Uy\h\^dyds 
Jo Jo Jo Jo 

fo fo ^^^!^^^^^ + ^^^) /* I I'^l I' -([0,1]) (/' 1^1'^^) 

+ e f [ \h\^dyds + C{e) [ [ uldyds + C{T) 
Jo Jo Jo Jo 

<6 [ [ ^-^^^^^dyds + € [ [ m^'dyds + CiT). 
Jo Jo ^ Jo Jo 

Substituting fl3.15p into the above inequality, we deduce (13.131) and (13.141) by 
taking sufficiently small e. □ 

In the following, we first verify that the temperature is dominated by 
velocity and magnetic fields. 

Lemma 3.9. One has 

e^dy+ [ [ nO^eldyds 
Jo Jo 

< C{T) + C{T) ( f f w^dyds) ' + C{T) ( f f h^dyds 
\Jo Jo / \Jo Jo / 

Proof. Multiplying 6^ on both sides of (I2.2ep . and then integrating it over 
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[0, 1] X [0,t], it leads to 

[ 9^dy+ I [ KO^eldyds 
Jo Jo Jo 



< 



< 



^ lo lo ^ ^^^y'^^y^^ ^ / / ^ ' ^^^^^^ 

C + C [ [ {d^^ + \uy\^)dyds + C f {\Wy\\l^ + \\hy\\lJ) f e^dyds 
Jo Jo Jo Jo 



t /.I 



< C + C 




{6^^ + \uyY)dyds 



+ C 



Jo 

t rl 




Wyydyds ] + C 



JO 



t fi 




JO 



hlydyds 



On the other hand, we get from fl2.2bp 



ht = ^hyy-{p+\\h\''), 

h\t=Q = ho{y), 

h\y=o^i = 0, 



(3.16) 



where h = ud^. The standard U estimates of solutions to linear parabolic 
problem yield 



t rl 




t rl 



\Uy\^dyds 



JO 




\hyy\^dyds 



'0 JO 

<c(i- 



"yyi 
-t ri 




JO 



< C{T) 1 



+ \h\')dyds 

t 




e^^dyds 



JO 
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and then 

•1 



t rl 








^0 



<c + c 
+ c 




2 



^yydyds 




^0 



e^'^dyds + C 



C 



ff 

'0 ^0 



hlydyds 



t /.I 



^0 




Uyl^dyds 



^0 



< C(T) 1 



t fi 




^0 



e^^dyds ] + C 
• 1 



t /.I 




^0 



Wyydyds ] + C 



ff 

'0 ^0 



hlydyds 



< C{T) + C / maxr / e'^dyds 



[0,1] 








Wyydyds ) + C 



^0 



t rl 




hlydyds 



^0 



It finishes the proof by Gronwall inequality and (13.8 



□ 



The following lemma also declare a relationship of density and velocity, 
as well as magnetic fields. 

Lemma 3.10. For any e > 0, it satisfies that 



t /.I 



vldy 




Ovldyds 



JO 



< C{T) + e 



t /.I 




wlydyds ] + £ 



JO 



t /.I 




^0 



hlydyds 



Proof. We can rewrite the equation fl2.2bp as follows 



u 



(3.17) 



Multiplying (m — -Vy) on both sides of (13.1 7p and integrating it over (0, 1) x 
(0, t), we find that 



t /.I 




ROuv. 



V 7 7 

-dyds 



(3.18) 



JO 




Jo 



- + 1 ^« + b • b. 



M — —Vy ) (i?/(is. 
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To complete the proof, it only evaluates all terms on the right hand side of 
fICT]) . Particularly, 

ROuv,, 




'0 JO 

< e 

< e 

< e 



t rl 



-dyds 




'0 ^0 



Ovydyds + Ce / max6' 



[0,1] 



u^dy I ds 




10 JO 

ft rl 



(3.19) 



Ovydyds + C^ / max 6'(?/, s)(is 



/ / evldyds + C{T), 
Jo Jo 



for any e > and further 



'0 ^0 



V 3 



u Vy ) dyds 



< c 




t rl 



y . 



+ 



^0 




dyds + £ I I nO^Oydyds 




^0 



K 
t rl 



~'0 



u Vy ] dyds 

V 



<C + £ 




^0 



KO^Oydyds + I max 



* e^ + e^ f'f A , , , 

u Vy ) dyds. 







[0,1] K 







V 



(3.20) 



In addition, by f l3.13p . we also obtain 

^ A 




^0 



h ■ hy ■ lu Vy 1 dyds 



t rl 



< c 




t rl 



|b ■ hy\^ dyds + C 



JO 




t rl 



< C{T) + C 




Jo 



^0 
2 



u Vy 1 dyds 



(3.21) 



u Vy ) dyds. 



Finally, substituting (13719]) - fl3:2T| ) into fl3J8|l . it leads to 



u Vy ) dy 







V 



t rl 




Ov^dyds 



Jo 



t rl 



< C{T) +e / nO^e'dyds 



+ 



C I ( max 

Jo ViO'i] 



Jo 
9^ + 9^ 



K 



+ 1 



u — —Vy ) dyds 
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which ends the proof by Gronwall inequahty and Lemma 13.91 . □ 

Combining Lemma 13.91 and Lemma I3.10[ we can give a priori estimates 
of velocity and magnetic fields, which play an important role to include the 
constant heat conductivity. 

Lemma 3.11. It holds that 

||b||L-((o,i)x(o,T)) + / + |wj^|^)dy 

.1 (3-22) 

+ / / {\hy\^+\Wy\''+\hyy\''+\Wyy\'')dydS<C{T). 

Jo Jo 

Proof. Multiplying f l2.2cp by w^^ and integrating it over [0, 1] x [0, t], one has 
\[ \'^y\'^dy = ^l \wy{y,0)\^dy- I j (h+^^^ ■ Wyydyds 

<C-C{T) / / \wyy\^dyds + C / {\hy\ + \vy\\wy\)\wyy\dyds 
Jo Jo Jo Jo 



n\wyy\'^dyds + C{T) / / \hy\'^dydi 
Jo Jo 



+ C(T) / max|w,,P( / vf,dy ] ds 
Jo [0.1] \Jo J 

<C-^^ f ['\Wyy\'dydS 



^ Jo JO ^ 

+ \ £ ( I "^lydyds^ ^ ^ (/ / Kydyds^ ) ^ (/ / "^yy^yds 



<C- ^ ' " ' / / \wyy\'^dyds + e / / hyydyds. 
^ Jo Jo Jo Jo 

Thus, 



/ |w,yp(iy+ / / \^yy'^dyds <C ^ e / / hyydyds, (3.23) 
Jo Jo Jo Jo Jo 

On the other hand, we can also rewrite (]2.2dp as follows 
ht = + - w H 

V V \ ^ ' y 
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1 

2 



Following the same procedure as above, we get from (13.41) 
\hy\'^dy 

<C-C [ [ \hyy\^dyds 
Jo Jo 

+ C / (|Mj^||b| + |wj/| + |t;y|m|)|byj^|(i?/(is 
Jo Jo 



<C — f I \hyy\^dyds + C sup |bp ( f Uydyds 

4 Jo Jo (i;,t)e{0,l)x(0,t) Jo Jo 

+ C max|b„P( / vf,dv]ds 

Jo m \Jo ^ / 

<C-—j J \hyy\^dyds + C j \h\^dy + - j \hy\'^dy 

<C{T)-^ [ [ \hyy\^dyds + ]- [ \hy\^dy + e [ [ wl^dyds, 
Jo Jo 4 Jo Jo Jo 

which implies 

/ \hy\'^dy+ / \hyy\^dyds<C{T)+e j I wl^dyds. 
Jo Jo Jo Jo Jo 



This together with (13.231) leads to 



[\\^^\^ + \hy\')dy + f [ 
Jo Jo Jo 



t rl 

i\Wyy\^ + \hyy\^)dydS < C{T). 



Similarly, we get 



4j„.j„ ^ / |2 / / |„, |2 




Wy\ dyds < I max|wy| | j \wy\ dy J ds 



'0 Jo Jo 

< 



C{T) [ I {\Wy\'+\V,yy\^)dydS 

Jo Jo 
< C{T). 
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and 

ft ^1 rt / A-l 

|2 



^0 



hy\ dyds < I max|bj/| ( / \hy\ dy ] ds 



<C{T) [ [ {\hy\^+\hyy\^)dydS 

Jo Jo 

< C{T). 

This completes the proof of the lemma. □ 

With fl3.22p in hand, we thereby can check Lemma 13.91 and Lemma 13.101 
again and deduce that 

Corollary 3.12. We have 

[ vldy+ [ [ 9vldyds<C{T). (3.24) 
Jo Jo Jo 

and 

I max0'?+^3rfs+ I 9^dy+ [ [ nO^Oldyds < C{T). (3.25) 
Jo [O'l] Jo Jo Jo 

Similarly, using estimates of solutions to linear parabolic problem 

again, we get from (13.141) and (13.161) 



t pi rt i-l 

Uy\^dyds = / \wyy\'^dyds 
^0 Jo Jo 

-t rl 

<C( 1 + 



f [ {p^ + \h\^)dyds 
Jo Jo 



t rl 



<c^i + J J e'^ds^ ^^-^^^ 

<c(l + [ max 9^ [ 9^dyds 
V Jo [0.1] Jo 

< C{T). 

In order to get the upper bound of temperature, we also need to establish 
higher order a priori estimates of {v,u,6,w,h). Thus, we introduce some 
auxiliary new variables motivated by |13] . 

X := [ [ {1 + 9'')9^dyds, (3.27) 
Jo Jo 

Y := max [ (1 + 9'^'')9ldy, (3.28) 

0<t<T Jg / 2/ ^' 



1 

2 

yy 



max / u,,,,dy. (3.29) 

0<t<T Jo 
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By interpolation and embedding theorem, It follows from ( I3.29p that 

<C(1 + Zi), (3.30) 

where | ■ | "-^^ = sup | ■ | . 

Lemma 3.13. Under the assumptions of Theorem \2.1\ we have 

max eiy.s) <C + CY^i^), (3.31) 

[0,l]x[0,t] 

for {y,t) G (0,1) x (0,T). 

Proof. By the embedding theorem, we deduce that 



max9''+\y,t) <C [ O'^+^dy + C [ (1 + 9^+% 

m Jo Jo 



y\dy, 



which implies by (13.251) and Holder inequality. 



[0,1] ^ 



<Cmax0^+i/" 9^dy + C [ {1 + 9y\9y\{l + O^dy 
[0,1] Jo Jo 

< ^max^«+^ + C + efWldy"^ ' + Ofdy"^ ' + C 

< I max ^^+5 + C{T) ( [ {I + eWdy] ' + C{T) 
^ [0.1] \Jo J 

< - max 9'^+^ + Cy 5 + C(T) . 
~ 2 [0,1] ^ ^ 

It finishes the proof. □ 

Lemma 3.14. Furthermore, we have 

X + Y <C{T) (l + Z^j . (3.32) 

Proof. We introduce the function as in 



Jo V 
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The simple calculation leads to 

Kt = + K,uy, (3.33) 

Kyt = (^~^y)^ + KyyVyUy + K^Uyy + ^^^^ Vy9t, (3.34) 

\K,\,\K.v\<C{l + 9'^+'). (3.35) 
Multiplying fl2.2ep by Kt and integrating it over (0, 1) x (0,t), we get 

Jo Jo \ V V V J 

+ 11 -OyKytdyds = 0, 
Jo Jo 

or equivalently 




JO 



dyds + / —6y ( —Oy] dyds 

Jn Jn V \ V 



V Jq Jq V ' \v "Jt ^ (3.36) 



= h + l2 + h + h + h + h, 

where /j, (i = 1, ■ ■ ■ , 6) is defined as follows 

Ji = - / / eeOtK^Uydyds 
Jo Jo 




^0 



V V V I V 



t /•! / \„2 12 



2/ 





'0 JO 

-OyKyyVyUydydS 



K^Uydyds 



'0 Jo 

h = - J j '^OyKyUyydyds 



n-6y (-) VyOtdyds. 



In the sequel, we will evaluate all terms f l3.36p according to fl3.33p - fl3.35l) . 
Firstly, one has by and fl3:28|) 

t rl a2 rt j-l 

'-dyds >C {l + e^){l + e'')eldyds > CX, (3.37) 



22 




^0 ^0 ^0 



and 



t fi 




^0 



—9.11 ( — 6'^ ) dyds 
V \v J t 



K 



7.^) dy-7; (-Oy) iy,m 



K 



(3.38) 



> c ^ (1 + e^feldy -c>CY-c. 



Secondly, we have by Cauchy-Schawtz inequality and (13.311) 

egOtKyUydyds 




^0 
t rl 



< c 




9+4 



Uy\dyds 



JO 



< 



eX + Ce [ / (1 + ey+^uldyds 
Jo Jo 



(3.39) 



<eX + C,YW) 



< six + Y) + C, 



t /.I 




Jo 



Uydyds 



for any fixed sufficiently small e > 0, and similarly, by recalling (I3.22p and 
(13. 26 p . one has 



V V 



-Otdyds 



1 1 

£ C- / / [(1 + + (1 + + Iw,!' + |b„p)] dydi 

Jo Jo 

<eX + C,\{l + 9Y+T^ [ [ uldyds 

Jo Jo 



(3.40) 



+ C,|(1 + 5)^|(°) f f 
JO Jo 

< e{X + Y) + C. 



+ \^y\'^ + \^y\^) dyds 
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and 



(Opguy 
^0 V ^ 
t fl 




\Uy /i|wyp z/|b 



y\ 

V 



' K^Uydyds 



< 



c f f {1 + ey+^ldyds + [ [ {i + ey+^\uyfdyds 

Jo Jo Jo Jo 

+ C [ [ il+ey+'\Uy\i\Wy\^ +\hy\^)dydS 

Jo Jo 

< + / / uldyds 

Jo Jo 

+ |(l + ^).+i|(0) f f lu^i^dyds 
Jo Jo 

+ C [ [ (1 + ey+'\Uy\{\Wy\^ + \hy\'')dydS 

Jo Jo 

ft fl 

<CYW1+CY^ / {ul + \wy\^ + \hy\^)dyds + C 

Jo Jo 

<eY + C. 

Thirdly, it also yields in view of (13. 8p and f l3.30p 




— 9y VyUydyds 

JO 

<C{ ^dyds) \{f [\l + ef^^\lvlc 




JO 



3(9+5) 

< eY + CZ^(^, 



JO 



and 



(3.41) 



(3.42) 



1/5 



t /.I 




K 



—dyK^Uyydyds 

10 JO 



<C( I I -^dyds 




JO 



< eY + CZ^. 



t fl 




JO 



(l + ef^+^^ulydyds) (3.43) 
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Finally, one has 



/ / -^y [-] VyOtdyds 

Jo Jo V ^V/v 

{1 + eyivyOtldyds 

<£X + C,|(1 + ^)^|(°) Tmaxf^ 

Jo [0,1] V ^ 



< C 



t r-l 



JO 



vldy ds 



<eX + C^Y^(5+^ 



<eX + CeY^(^ 



r ( 

I max 

Jo [0,1] V 

!' ( 

I max 

Jo [0,1] V 



V 



v^dy ] ds 



ds. 



By the embedding theorem again, one has 

2 



f max — - as 
Jo [0,1] V ^ / 



t rl 



< c 




JO 



dyds + C 



t /.I 




^0 



K,9,„ ( kQ 



V \ V 



dyds 



^dyd 



+ C[ I I -^dyds 




JO 



' "'^O^f^) dyds 




^0 



<C\\{l + 9) 



g+2|(0) 



+ 



/ / {i + er'{eX + eV,ul + ul + \wX + \hy\')dyds 

Jo JO 

In particular, we also have 

r [\i + ey+'eXdyds 

Jo Jo 



< (71(1 + ^)81(0) f f ^i + oyo^dyds 
Jo Jo 
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and 

Jo Jo 

<c [ /" (1 + ey+^^ldyds 

Jo Jo 

<c\{i + eY+^T^ f [\ldyds 



JO 



g+10 

< C 1 + 



and 



r [\l + er+' {ul + \wX + M) dyds 
Jo Jo 

<C|(l + 0r+2|(°) / / {ul + \wy\^ + \h/)dyds 
Jo Jo 



< C ( 1 + 



which imphes 



y2(/+5) / max ( — 1 ds 
Jo [0,1] V ^ 



/ £+1 1 q+i 3g+10 ' 

< C 1 + F^+S + X2F2(9+5) _|_ y 4(9+5) 

<e{X + Y), 

which together with fl3.44p shows that 

\Ie\<eiX + Y). (3.45) 

This together with fl3.36p - fl3.43p completes the proof. □ 

With the above a priori estimates, we can succeed in obtaining the upper 
bound of temperature, which is another essential contribution of the paper. 

Lemma 3.15. The following a priori estimates hold, 

+ + <C{T). (3.46) 



and 



!\dl + + ul)dy + f f\dl + |bi|2 + ul)dyds < C{T), (3.47) 
Jo Jo Jo 
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Proof. Differentiate (I2.2bp with respect to t, multiply it by Ut, and then 
integrate to obtain 

As a consequence, it implies that 

ufdy + / v^y^dyds 
Jo Jo 



ul\uyt\dyds + C / \ptUyt\dyds + C / |b ■ btMyt|(i?/(is 



< 

< - [ [ ulfdyds + C [ [ utdyds 
~ ^Jo Jo " Jo Jo " 

+ C [ [ pldyds + C [ [ \ht\^dyds + C{T) 
Jo Jo Jo Jo 

ft rl rt rl rt rl 



<C(T) + ^ / / uf^tdyds + C I I pfdyds + C I / \ht\''dyds 



^2 ^„,j„ , ^ / / „2j„,j„ , ^ / / lU |2, 
'0 JO JO JO JO JO 



which leads to 

• 1 rt /•! 



Jo Jo 



ufdy + I I v^y^dyds 



<C{T) + C [ [ pldyds + C f [ \ht\^dyds. 
Jo Jo Jo Jo 



(3.48) 



We notice that 

rl rt rl 



npldyds = / / {p^vt + peOtf dyds 
Jo Jo 



t rl / j^Q / 

—Uy + f i?p + —e^ ) Qt ) dyds 

t rl rt rl 




Jo 



lo lo ^'''y'^y'^' ^ ^ / / ^ e^tdyds (3 49) 

<^(T)|^2|(o) /"* f\^^dyds + C f f\l + e^)e^,dyds 
Jo Jo Jo Jo 



< C{T)Y~' + X + XY~^ 

/ q+8 \ 

< C(T) (l + Z^ \ . 
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and follows from and flX^ 

\ht\^dyds 




^0 



< [ [ {\h\\l + \Wy\^ + \hyy\^ + \hy\%l) dydS 

Jo Jo 

<C I I uldyds + C I I {\wy\^ + \hyy\^)dyds 
Jo Jo Jo Jo 

+ C max|b„p| / vfdy ] ds 
Jo [0,1] ' ^' Vio ^ / 

C + C [ [ {\hy\'+\hyy\^)dydS 

Jo Jo 



< 

< C{T). 

Substituting fl3:49|l and fl33nD into ([33HD, it shows that 

1 ft fi 

2j„, I / / „,2 



/ g+8 \ 

< c (i + z^A . 



On the other hand, it is also deduced from fl2.2b|) that 



Uyy = - \Ut+ \ P + 



XVyUy 



\ ^ ^ y 

and then by integration, it follows from f l3.5ip 



(3.50) 



Utdy+ / U fdyds 

'o Jo (3.51) 



/ ulydy <C [ (m? +pI + |b|2 . \hy\^ + vlul) dy + C 
Jo Jo 

<C [ uldy + C f pldy + C [ v^uldy + C [ \hy\^dy + C 
Jo Jo Jo Jo 

<c(i + z^o+ f\i + eyidy + {\eT^ + \ul\^'^) fvldy 

\ Jo Jo / 

f 9+8 q+8 ; 

< C M + ZiTTo + z- 

< C{T) (l + 



9+8 
; q + W 



Thus we have Z < C due to < -2±^ < 1, and X and Y are also bounded. 

— <j+10 ' 

Subsequently, \9\^^\ /q (uf + Oy + Uyy)dy and /g {u1^ + 9^)dyds are 

also bounded. □ 
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In the sequel, the lower bound of temperature is obtained. 
Lemma 3.16. One has 

e{y,t)>CiT). 
Proof. Let = |, then fl3.12p is written as 



K 



2Kei /i|wj2e2 u\hy\^e^ 



_ ( _ , YEi 

V V' 2X&J + 4A' 

which implies 

e,<-(-0,) +c(T), 

by (13.21) and fl3.46p for some positive constant. Define the operator : = 
-i + -# \ and then 

< 0, on Qt = (0, 1) X (0, T), 

e\t=o>0 on [0,1], 

Qy\y=o,i = on [0,T], 

where Q{y,t) = C{T)t + max[o,i] — Q{y,t), and by the comparison the- 
orem, one has 

min Qjy, t) > 0, 

{y,t)eQT 



which is inferred 



1 



[0,1] ^o(y), 

for any {y,t) E Qt- □ 

3.3 Higher order derivatives a priori estimates of (w, b, 9) 

Lemma 3.17. One has 



(W„b„t;,)|(°)+ / (\^^,\^+\(vh)t\^+\Wyy\^+\hyy\^)dy 

Jo 

+ 11 {\^ty? + \hty\'' + ely)dyds<C{T). 
Jo Jo 



(3.52) 
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Proof. Firstly, we differentiate ( I2.2cp with respect to t, multiply it with 
and then integrate over (0, 1) x (0,t) 



1 , r 




JO 



V 



= \ I Wt?{y^^)dy+ [ [ ^UyWy ■ wtydyds - I I hfWtydyds 
^ Jo Jo Jo Jo Jo 

<C+\I [ ^\wty\^dyds + C{T) [ [ {ul + \wy\''+\ht\^)dyds, 
Jo Jo Jo Jo 

which follows from ^M), and iKAT} 

[ \wt\^dy+ [ [ \wty\^dyds<C{T). (3.53) 
Jo Jo Jo 

Similarly, we get from fl2.2cp 



which leads to 



/ \wyy\^dy<C [ {\wt\^ + Ihyl"^ + vl\wy\'^)dy 
Jo Jo 



<C + Cmax|w„| / v,,dv 

[0,1] ' Jo ' 

<C{T) + C f \wy\^dy + l- [ \wyy\'dy 
Jo ^ Jo 

<CiT) + ^l^ \wyy\'dy. 



I.e., 



1 

l^yyl^dy < C{T). 



Secondly, we further deduce that 



t /.I ^2 

, -i^ly^y^^ 

Jo ^ 




<C [ f {dl + ul + ul + |w,|^ + |b,|^ + vpl + el)dyds 



Jo 



<C + C [ maxOl [ {vl + dl)dyds 
Jo [0.1] ' Jo ' ' 

<C + C / eldyds +- -Adyds, 
Jo Jo Jo Jo ^ 
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and then 

[ [ elydyds<C{T) + C f [ eldyds<C{T). (3.54) 
Jo Jo Jo Jo 

In addition, by f l3.5p . one has 

X{\nv)y =X{\nvo{y))y+ / {py + h ■ hy)ds + {u - uq) 

Jo 

or equivalently 

v'y<C + C [\\hy\' + + pleads 

Jo 

<c + c f [\\hy\' + \hyy\' + el + e'jdyds + c fvlds 

Jo Jo Jo 



t 

2. 



<C + C v'yds. 
Jo 

Thus, we deduce that Vy is bounded by Gronwall inequahty. Lastly, differ- 
entiate fl2.2dp with respect to t and then muhiply by {vh)t and integrate to 
get 

1 d 



I ivh)ldy+ / ulhytl'^dy 
Jo Jo 

- I -hyt ■ {hUyy + hyUy + y) dy 

Jo 

+ / —Uyhy ■ {hUyy + hyUy + bjUy + htyV)dy 

Jo 

- Wf (huyy + hyUy + ht'Vy + htyv)dy 
Jo 



(3.55) 











Jo 





So we get from flX^ . fIXiTD and flX^ 

\ {vhftdy+ [ [ \hyt\^dyds<C, (3.56) 
Jo Jo Jo 

by integration for any fixed sufficiently small e. Furthermore, by (12.2dp . one 
has 

\hyy\'dy < C f{{vh)] + |w,|2 + vl\hy\^)dy 
Jo Jo 

<C{T) + C [ \hy\^dy<C{T), 
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and similarly 



\hy\^<cf \hy\^dy + C f \hyy\^dy<C{T). 
Jo Jo 



This completes the proof. 
Lemma 3.18. We have 



\Oyf^+ [\o^ + Oly)dy+ f f eldyds<C{T). 
Jo Jo Jo 



□ 



(3.57) 



Proof. Differentiate fl2.2el) with respect to t, multiply it by eo9t and then 
integrate it over [0, 1] 



dtJo 2 

• 1 



dy+ I -eeOltdy 



- peeeUyOt - dpg^ceu dt - OpeeCeUyOt - dpgCeUytOt 



+ —egUyUytOt - ^egUyOt + —ee^y ■ ^ym - -^egUylWyl Ot 
H eeOthy ■ hyt 



V 



V/ V 



dvVyUydyOt 



K 



V / V 



eeeUyOyOt 



K 



K 



egUyOyO. 



yt 



eovVyOyOl 



2q2 



— ( ^^edtdydyt — —eevVy9t9yt — —eeeOyOtOyt 



eeedydt 



which implies by integration 



t /.I 



efdy + 

< C (1 




'0 Jo 



Oy^dyds 



[ [ {vl + el)eldyds 

Jo Jo 



<C{1 



max O'^ds 

[0,1] 



t /.I 



<C[l+e 




t 



Oyfdyds + 



JO 




O^dyds 



JO 



and then O^dy and /* Oy^dyds are bounded. Taking the same operation 
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as hyy and w^^, we deduce from (I2.2ep 
Jo Jo 1^ 

<e f elydy + C{T), 
Jo 

which imphes the boundedness of jj^ ^yy^U} ^"^^ so is |6'j,|*-°^ by employing the 
embedding theorem. □ 

Now we have estabhshed Lemmas 13.11 - 13.18] we can prove that the Holder 
estimates of solutions by the routine manner. Indeed, from (13.461) . f l3.52p and 
fl3.57p . one has 

Kw,b,^)GCi'°(gT)'. 

Following the procedure of [2], [22] or [28], we can get 

(^i,w,b,^)GCi'^(gT)', 

Finally, by the classical schauder estimates and the methods in [TJ], the 
Holder estimates of solutions are derived. This completes the proof of The- 
orem [2]TJ 
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